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Resolving the coronal line region of NGC 1068 with near-infrared 
integral field spectroscopy* 
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ABSTRACT 

We present adaptive optics-assisted J- and /iT-band integral field spectroscopy of the inner 
300 X 300 pc of the Seyfert 2 galaxy NGC 1068. The data were obtained with the Gemini 
NIFS integral field unit spectrometer, which provided us with high-spatial and high-spectral 
resolution sampling. The wavelength range covered by the observations allowed us to study 
the [Caviii], [Sivi], [Si vii], [Alix] and [S ix] coronal line (CL) emission, covering ionization 
potentials up to 328 eV. The observations reveal very rich and complex structures, both in 
terms of velocity fields and emission-line ratios. The CL emission is elongated along the NE- 
S W direction, with the stronger emission preferentially localized to the NE of the nucleus. CLs 
are emitted by gas covering a wide range of velocities, with maximum blueshifts/redshifts of 

1600/1000 km s"'. There is a trend for the gas located on the NE side of the nucleus to 

be blueshifted while the gas located towards the SW is redshifted. The morphology and the 
kinematics of the near-infrared CLs are in very good agreement with the ones displayed by 
low-ionization lines and optical CLs, suggesting a common origin. The line flux distributions, 
velocity maps, ionization structure (traced by the [Si vii]/[Si vi] emission-line ratio) and low- 
ionization emission-line ratios (i.e. [Fe n]/Pa/? and [Fen]/[Pn]) suggest that the radio jet plays 
an important role in the structure of the CL region of this object, and possibly in its kinematics. 

Key words: galaxies: active - galaxies: nuclei - galaxies: kinematics and dynamics - line: 
formation - line: profiles - infrared: galaxies 



1 INTRODUCTION 

Feedback due to outflows and jets is believed to be a major driver 
of the co-e volution of ga laxies and their supermassive black holes 
(review bv lFabianll2012h . Radio jets interact in various ways with 
their surroundings, sweeping up, and driving shocks into the ambi- 
ent gas. While direct jet-gas interaction locally affects the velocity 
field of the surrounding medium, shocks or the radiation field from 
fast, auto-ionizing shocks can significantly influence or dominate 
emission-line strengths. Nearby galaxies are key laboratories for 
studying jet-gas interaction in great detail, and at high resolution. 
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through spatially resolved spectroscopy of the (high-ionization) 
emission lines. 

A subset of active galactic nuclei (AGN) show in their spectra 
emission lines from very highly ionized atoms, known as 'coronal 
lines' (CLs). They are coUisionally excited forbidden transitions 
within low-lying levels of ionized species with ionization poten- 
tials IP > 100 eV. CLs have been detected in the optical and in- 
frared spectra of all types of AGN, including Seyfert 1 and Seyfert 
2 galaxies, narrow-line Seyfert 1 galaxies, and radio galaxies 



(e.g. Penston et al. 


1984; Marconi et al.' 1994; Nagao et al. 2000|; 


Sturm et al. 2002; 


Rodriguez- Ardila et al.. .2002. 2006; ,Deo et ak 


2007: MuUanev & Wardll2008l; Komossa et alj 2008l;lGelbord et alj 


2009i). Their abundance is aroroximatelv similar in AGN of type 1 



and 2 (Rodriguez- Ardila et al. 201_1) and they represent one of the 
major gaseous components of active nuclei. 

The precise nature and origin of CLs are still a matter of 
debate. Different scenarios have been considered i n the litera- 
ture, including winds from the molecular torus (e.g. IPier & Voi3 
ll995l ;iNa gao et al. 2 00 0; MuUanev et al. 20091 (X-ray) ionized ab - 
sorbers (NLR; e g iKomossa & Finklll997al lbl; |Porauet et alj|l999h . 
a high-ionization component of the inner narrow-line region 
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Komossa & Schulj[T997l : iFerguson efZII 19971 : Isinette et al] 



1 1997), and a low-density component of the interstellar medium 



it y 

(ISM: iKorista & FerlandlT989i) . However, up to now none of the 
above hypotheses has been fully discarded because of the lack of 
high-spatial resolution observations on a sizeable number of tar- 
gets. 

The mechanisms powering the CLs have also been a subject of 
strong debate, with two competing processes usually invoked: pho- 
toionization from the central engine and shocks be tween the jets 
and th e NLR gas. Dopit a & Sutherland! ( Il996h and lBicknell et alj 
( Il998h proposed that the emission in the NLR may be entirely 
caused b y shocks. Velocity splitting of over 1000 km s"', re- 
ported by lAxon et al.l ( Il99 8'). in the vicinity of some of the bright 
emission-line knots in the NLR of NGC 1068 provides evidence of 
the existence of fas t shocks in AGNs. However, Hubble Space Tele- 
scope (HST) da ta dCrenshaw & Kraemej|200ol : ICecil et al.ll20o3 : 
iMazzalav et"ai]l2010l) have shown emission-line ratios consistent 
with photoionization instead of shock heating mechanisms. More- 
over, models accounting for both photoionization from the cen- 
tral radiation source and shocks were required to explain both 
the continuum and NLR spectrum, includin g the CLs, observed 
in the spectra of a sample active galaxies ( Contini et al.l |2003I : 
iRodrfguez-Ardila et al.l 12006). Furthermore, sometimes the same 
data sets have been explained by different scenarios. In particu- 
lar, both mechanisms, shocks and photoionization, may plausibly 
be at work in the same galaxy. Therefore, high-resolution studies 
of nearby galaxies are essential in order to disentangle the effects 
and contributions of shocks and photoionization. 

Recent advances in the study of CLs include the determina- 
tion of the size and morphology of the region emitting these lines 
by means of optical HST and ground-based near-infrared (NIR) 
adaptive optics (AO) ima ging/spect roscopy in a few AGNs (e.g. 



Riffel 



ve optics (P 
et alj l2008l: 



Storchi-B ergmann et i 



op y m a 
ai]l2009l: 



20 id : iMuller-Sanchez et al.,201 1.) . iMazzalav et al. 



Mazzalay et alj 



j2010h , for in- 
stance, presented an analysis of the optical CLs of a sample of 10 
Seyfert galaxies, including NGC 1068. These authors used the ex- 
cellent spatial resolution provided by STIS/HST to resolve the CL 
region (CLR) and study the properties of the CLs as a function of 
the distance to the nucleus. The optical CLs were found to be very 
complex, showing strong variations in their intensities and profiles 
along the spatial direction, which in most cases was nearly aligned 
to the position of the radio jet. These lines are emitted from very 
compact regions (~ 30 pc) up to regions of ~ 250 pc. The mor- 
phology and the kinematics of the lower-ionization CLs (i.e. [Ne v], 
[Fe vii] and [Fe x]) are very similar to that of lower ionization lines 
such [Oiii]. On the other hand, CLs with IP > 250 eV (i.e. [Fexi] 
and [S xii]) are emitted in much more compact regions and do not 
seem to share the same overall kinematics found for the lower ion- 
ization lines. Several pieces of evidence (e.g. lack of correlation 
between the radio and CL emission, and relatively good success of 
photoionization models in reproducing measured line ratios) point 
towards photoionization from a central source as the main mecha- 
nism of formation of CLs. 

However, the STIS results are restricted by the relatively small 
region covered by the slit width and along a particular direction that 
can be analysed with long-slit spectroscopy. In spite of the typical 
0.1 arcsec spatial resolution, the limited area covered by the slit 
did not allo w rejection of jet-NL R gas interaction in the production 
of the CLs ( IMazzalav et alj|201 0). A further step on this subject is 
expected to be achieved by means of AO integral field spectroscopy, 
which can give us a more detailed view of the innermost parsecs of 
nearby AGNs at similar angular resolutions as HST. 



NGC 1068 is a natural candidate for a follow up with AO inte- 
gral field spectroscopy. This object is one of the nearest and proba- 
bly the most intensely studied Seyfert 2 galaxy in the literature. Ob- 
servations in all wavelength bands from radio to hard X-rays have 
formed a uniquely detailed picture of that source. It hosts a promi- 
nent NLR that is approximately co-spatial w ith a linear radio source 
with two lobes jWilson & Ulvestadll983h , extensively character- 
ized froni_subarcsecord_c^ probed by the HST jEvans et al.l 
I1991I : iMacchetto et al. || l994|b to the ionization cone and extended 
emission-line region ( !Poggei l98 8l: IUngeretal]|l992h extending to 
radii of at least 30 arcsec (1 arcsec = 73 pc for z = 0.003793 and 
H„ = 75kms-' Mpc"'). 

The op tical and NIR CLs of NGC 1068 were first studied by 
lMarconi"et al. ( 1996). They found that the flux distribution peaks 
50 pc NE of the nucleus and appears to arise predominantly in 
outflowing gas within the prominent ionization cone aligned with 
the radio jet. The line ratios argue against collisional ionization but 
could be explained with photoionization by radiation from the ac- 
tive nucleus. Prieto et al. (2005) presented a [Si vii] 2.48 yum narrow 
band image obtained with the IR camera NACO at the ESO/VLT 
(spatial resolution FWHM = 6.8 pc or 0.097 arcsec). The im- 
ages show a diffuse coronal region with filamentary structure. The 
extent of the CL emission is larger than that predicted by pho- 
toionization models, which argues for additional in situ gas ex- 
citation, the most plausible en ergy source being shock excitation. 
iRodrfguez-Ardila et al. I ( I2OO6I) reported a double-peak structure in 
the [Fe vii] /16087 line, with the relative separation between the two 
peaks increasing with the distance from the centre, signalling the 
presence of a highly energe tic outflow, not detected in l ower ion- 
ization lines. Very recently. iMiiller-Sanchez et"ak l (l20Tl hereafter 
MSll) analysed the morphology and kinematics of [Sivi] 1.96 yum 
(IP=166 eV) CL of NGC 1068 by means of /f-band AO integral 
field spectroscopy SINFONI data. These authors confirmed the 
great complexity displayed by the [Sivi] emission-line region of 
this galaxy and favoured a scenario where the radio jet has no strong 
influence in the formation and kinematics of the CLR, probably 
only affecting the kinematics of some emission-line knots located 
at distances of less than 1 arcsec from the nucleus, towards the NE. 

Because of the above points, additional information is needed 
to remove these ambiguities and identify the dominant mechanism 
powering the CLs. Moreover, the role of the radio jet in shaping 
the morphology and kinematics of the inner NLR gas can now be 
studied in greater detail thanks to the availability of spatial reso- 
lution at the parsec scale using ground-based facilities. Here, we 
present AO-assisted /- and K-b&nA integral field spectroscopy data 
of the 300x300 central parsecs of NGC 1068 obtained with NIFS at 
Gemini North. The wavelength range covered by these data (1.15- 
2.5 yum) allowed us to study in a spatially resolved fashion the 
[Caviii], [Sivi], [Sivn], [Alix] and [Six] CL emission, covering 
in this way ionization potentials up to 328 eV. In Section |2l we 
describe the observations and data reduction procedures. The anal- 
ysis of the morphology and kinematics of the CL gas is described 
in Sections [3.2l and [373] respectively. Section l34l describes the ion- 
ization structure of the CLR and the analysis of the origin of the 
[Fe 11] emission in this galaxy. The summary and final conclusions 
can be found in Section|4] 



2 OBSERVATIONS AND DATA REDUCTION 



Spect ra of NGC 1068 were obtained with NIFS dMcGregor et al.l 
I2OO3I) on the Gemini North telescope in December 2006 under pro- 
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gramme GN-2006B-C-9. T he observations and data reductio n have 
been described in detail by IStorchi-Bergmann et all (l2012h . NIFS 
was used with the ALTAIR facihty adaptive-optics system in its 
natural guide star mode. The compact Seyfert nucleus of NGC 1 068 
was used as the adaptive-optics reference object. The J, /fsijo,., and 
A'long observations discussed in this paper were obtained on 2006 
December 15, 12 and 13 (UT), respectively. The y-band spectra 
have a resolving power of ~ 6040, and the A^-band spectra have a 
resolving power of ~ 5290. 

NIFS is an image-slicer integral-field spectrograph with a 
square field of view of 3 x 3 arcsec^, which is divided into 29 slitlets 
each 0.103 arcsec wide with a sampling of 0.044 arcsec along each 
slitlet. The instrument was set to a position angle on the sky of 300° 
to align the slitlets approximately perpendicular to the axis of the 
radio jet and ionization cone observed in the nucleus of NGC 1068. 
This provided coarser spatial sampling along the jet axis, and finer 
spatial sampling across the jet. 

Each data set consisted of sets of multiple 90 s exposures 
that began with an offset sky position and were followed by nine 
galaxy exposures obtained in a 3 x 3-frame grid centred on the 
galaxy nucleus. A frame offset of 1.0 arcsec was applied along and 
peipendicular to the NIFS slitlets, so that the full field of view is 
5.0 X 5.0 arcsec- on the sky, with maximum exposure in the cen- 
tral 3.0 X 3.0 arcsec^. This exposure grid was repeated nine times 
for each grating setting. The nucleus of NGC 1068 was within each 
frame, so it was used to spatially register the exposures. An arc 
lamp exposure and a measurement of a telluric standard star were 
obtained after each data set. 

The data reduction was performed using the gemini iraf 
NIFS package and followed sta ndard procedures as described by 
IStorchi-Bergmann et al.l ( |2012|) . The data were flux calibrated us- 
ing the observations of the telluric standard star, HIP 5886, which 
was assumed to have the intrinsic spectrum with a level set by its 
2MASS magnitude in the wavelength band of the spectrum and a 
shape defined by a 8600 K blackbody, which has the same 2MASS 
J — K colour as HIP 5886. The data reduction resulted in single 
data cubes at each grating setting with all instrumental signatures 
removed. 



3 RESULTS AND DISCUSSION 
3.1 The CL spectrum of NGC 1068 

The NIR spectrum of NGC 1068 displays numerous emission 
lines, the brightest of which are Pa^S, [Feii] 1.26 /im, several H2 
lines, and the high-ionization CLs [Six] 1.25 yum (IP=328 eV), 
[Sivn] 2.48 yum (IP=205 eV) and [Sivi] 1.96 urn (IP=166 eV). 
Additionally, NGC 1068 exhibits weak CL emission from [Alix] 
at 2.04 yum (IP=285 eV) and [Cavm] at 2.32 /im (IP=128 eV). 
Fig.[T]shows the 7, A'shon and A'long continuum-subtracted spectrum 
of NGC 1068 integrated in a region size of ~ 0.2 x 0.2 arcsec- 
located at a distance of ~ 0.6 arcsec north from the nucleus. It 
can be seen that the CLs virtually dominate the NIR spectrum. The 
strongest CLs observed, [S ix], [Si vi] and [Si vii], are blended with 
the [Fell] 1.26 yum, H2 1.96 yum and Ht 2.48 /jm emission lines, 
respectively. Therefore, in order to study the CL emission it was 
first necessary to subtract the contribution from the [Fe 11] and H2 
lines. For this purpose, we first used the continuum task of iraf to 
subtract the continuum emission from the spectra of the data cubes. 
The simplest case to handle was the [S ix]-H[Fe 11] blend, since these 
two lines are not very strongly mixed. For each spectrum, we rep- 
resented the profile of each of these lines with a Gaussian function. 



and subtracted the [Feii] emission from the data. This procedure 
was satisfactory for most of the spectra, although in some regions 
it was necessary to include an extra Gaussian component to ac- 
count for the strong asymmetries or double peaks displayed by the 
[Fell] line. These regions are located towards the SW of the nu- 
cleus, where the [Fe 11] line splits into two components, generally 
one redshifted and one blueshifted (up to ~ -450 km s"'). Note 
that, even though the [Feii] show large blueshifts, the separation in 
wavelength with the [S ix] line is still large enough to allow a clear 
deblending of the lines. Detailed analysis of the kinematics of the 
[Fe 11] emitting gas, based on this same data set, can be found in Rif- 
fel et al. (in preparation). Deblending the [Si vi]-l-H2 lines was more 
challenging than the previous case because of the strong broad and 
complex profile displayed by [Si vi] in many regions, making the 
discrimination of each line difficult (see Fig. [T}. For this task, we 
used as reference for H2 1.96 yum the H2 2.12 /im line. The latter 
was chosen because it is located relatively close in wavelength to 
the former and in a region free of any other features. We modelled 
the H2 2.12 yum line as the sum of two Gaussian components (to 
account for the asymmetries or double peaks displayed by this line 
in some regions), and assigned the same widths, velocity shifts and 
relative peak intensities of these two components to the ones used 
to describe the emission of H2 1.96 fim. Additional constraints were 
necessary in some critical regions where the width and the inten- 
sity of the [Si vi] line did not allow a clear discrimination between 
the coronal and molecular lines. In these regions, we applied an 
additional restriction, fixing the value of the H2 1.96 yum/2.12 yum 
line ratio to 1.08. This corresponds to the mean value derived from 
the measurements in the regions where the H2+[Sivi] blend was 
not too strong. Although this line ratio is not independent of the 
physical properties of the emission-line gas, the measured values 
showed no significant variations between the spectra. Once mod- 
elled satisfactorily, the H2 1.96 //m emission was subtracted from 
the spectra. 

As for the [Si vi] line above, the observed [Si vii] line is lo- 
cated in a spectral region containing H2 emission lines. Neverthe- 
less, because of the low intensity and particular distribution of the 
H2 emission (see Section [X2] l. its contribution to the [Sivii] emis- 
sion is practically negligible. Only when the [Si vii] is analysed in 
velocity bins (Section [3.3l >. a clear contamination is present in the 
region where the H2 emission reaches its maximum (~ 1 arcsec 
east of the nucleus). Since the H2 emission does not contribute sig- 
nificantly to the [Si vii] emission and, when it does, it can be easily 
pinpointed, we did not subtracted the H2 emission lines from the 
A'long NIFS spectra. 



3.2 Morphology and extension of the CLR 

After applying the cleaning process described above, we con- 
structed bi-dimensional maps of the flux distribution of the CLs ob- 
served in the /, A'shon, and Anions data cubes of NGC 1068. We also 
constructed flux distribution maps of other lines of interest, namely 
PayS, Bry, H2 2.12 yum, and [Feii] 1.26 yum, for comparison pur- 
poses. The maps were constructed by integrating the flux under the 
emission-line profiles observed in the continuum-subtracted spec- 
tra. We used wavelength bins wide enough to include the emission 
from the highest velocity line-emitting gas. Figs.|2]and[3]show the 
flux distribution maps of the CLs and the selected low-ionization 
and molecular lines of the inner ~ 5 x 5 arcsec^ of NGC 1068. The 
maps are centred at the position of the continuum peak of the corre- 
sponding and A'-band, which coincides with the location of the 
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Figure 1. NGC 1068 continuum-subtracted spectra in the rest frame of the J (upper panel), S'shoit (middle panel) and S'long (lower panel) bands extracted from 
a ~ 0.2 X 0.2 arcsec" region centred at ~ 0.6 arcsec north of the nucleus. The labels mark the position of the principal emission lines observed in the spectra. 



central black hole jThompson et alfcoOlh . In these maps, we only 
report values with a signal-to-noise ratio (SNR) higher than 2. 

Visual inspection of Fig. |2] shows that the region emitting the 
CLs of NGC 1068 is highly inhomogeneous and elongated along 
the NE-SW direction. The inhomogeneous distribution can also be 
seen in the flux distribution maps of the low-ionization lines (i.e. 
Pa/3, Bry and [Feii], Fig.[3}- The [Sivi] map presented here is in 
agreement with the one reported by MS 11. However, some dif- 
ferences are observed due to the fact that these authors did not 
subtract the H2 1.96 yum emission, which can severely contami- 
nate the [Si vi] map at some particular locations. Although the inte- 
grated flux distribution is not strongly affected by this, the analysis 
of the velocity channels is more compromised, especially in the 
high-velocity blueshifted bins (see Section [331 for details). The ex- 
tension of the [Six], [Sivi] and [Sivii] CLs towards the NE from 
the nucleus is very similar, extending to a distance of 2.3 arcsec 
(~ 170 pc). Note, however, that the NIFS data do not cover regions 
beyond this. Therefore, we cannot exclude the possibility of emis- 
sion at larger distances from the nucleus. [Ca viii] and [Al ix] are 
observed in a more compact region, extending up to ~ 140 pc NE 
from the nucleus. Towards the SW the CLs are weaker, displaying 
an abrupt decrease in their intensity. [Si vi] and [Si vii] show a sec- 
ondary peak of emission at ~ 2.2 arcsec SW, which can also be 



seen in the lower-ionization lines of Pa/?, Bry and [Feii] (Fig. |3]l- 
This feature is also present in the data presented by MS 11, con- 
firming that it is real and not an instrumental artefact. Although the 
extension of [S ix] towards the NE is very similar to that of [Si vi] 
and [Si vii], towards the SW this line is observed only up to a few 
arcseconds from the nucleus, and no emission is detected at the po- 
sition of the emission peak observed at ~ 2.2 arcsec SW in the other 
CLs. Note that, when comparing the spatial extension of different 
lines one must keep in mind the SNR at which these lines are de- 
tected, as the different sizes could be a simple artefact of the lower 
SNR of the weaker lines. 

At this point, it is useful to compare the CL emission found 
here with the continuum and line distributions at other spec- 
tral regions. NGC 1068 harbours a nuclear jet observed in ra- 
dio, which could be playing an important role in the ionization 
structure and kinematic s of the NLR and, in particular, of the 
CLs of this galaxy (e .g. iMacchetto et all 1 1994 1 Axon et al.lll998l; 
iMazzalav et alj[2010l) . Panels a-c of Fig. |4| show a superposition 
of the [Sivi] CL map and the 6 cm radio emission observed at 
high- and low-resolutions with VLA and MER LIN, respectively . 
The radio data were taken from the work of Gallimore et al.l(ll996t) . 
Moreover, it is possible to compare the infrared CL emission with 
that of the optical [O iii] 5007 A line at similar spatial resolutions 
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Figure 2. Flux distribution maps of the CLs observed in the / and K data cubes. Each line is indicated in the upper right comer of the con'esponding panel. 
The colour bars indicate the values of integrated flux in units of 10"'^ erg cm"^ sec"' . The orientation and spatial scale (indicated in the upper left panel) are 
the same for all the panels. The position of the continuum emission peak is indicated in each panel by a plus symbol. 



by means of the F502N WFPC2/HST image jPressel et al1ll997h . 
The [O m] image was taken from the HST archive and its contours 
are shown in panels d and e of Fig.|4] superimposed on the [Si vi] 
map. For the superposition of the different data sets we assumed 
that the maximum of the NIR continuum and the Si knot of the jet 
observed with MERLIN coincide with the location of the central 
black hole, while the [Qui] B cloud is sl ightly offset towards the 
NE (see discussion in lGalliano et al.ll2003l and references therein). 

It is easy to see from Figs. |2] and |4] that the regions where the 
CLs are more intense are the ones located towards the N-NE from 
the nucleus, forming an arc that starts in the nucleus and extends 
about 0.8 arcsec towards the NE. This arc is resolved into individu- 
als clouds by the HST [O iii] image (clouds B, C, E and D of panel e 
of Fig.|4j. Noteworthy is the spatial coincidence between the [O iii] 
and CLs emission, both in large and small scales. The CLs peak at 
~ 1.7 arcsec NE can be associated with the cloud G of [O m], while 
cloud F also has a counterpart of CL emission except for [Alix], 
whose emission is not detected in this region. The only exception 
to this point-to-point correspondence is the secondary peak of emis- 
sion displayed by the [Si vi] and [Si vi] lines at 2.2 arcsec towards 
the SW, which is not observed in the [O m] emission. The lack of 
detection of this latter line in the optical is probably due to the 



strong extinction affecting NGC 1068 in the SW jMazzalav et alj 
[2010, and references therein). 

An inspection of the [S ix] map in Fig.|2]shows that the spatial 
distribution of this line is smoother and more symmetric than the 
other CLs. This can be explained by the poorer performance of the 
AO system at short wavelengths, which results in a loss in spatial 
resolution of the y-band data compared to ones at the A'-bands. The 
comparison of the flux distribution maps of Pa/3 and Bry (Fig. [3]l 
confirms this hypothesis: they are formed under the same physical 
conditions and whenever both are detected, they should exhibit a 
similar distribution (albeit a scaling factor). However, it is clear 
from Fig. |3]that Bry traces more clearly the NLR structure of this 
galaxy, while the Pa/3 distribution is more diffuse. 

The excellent match between the flux distribution of the NIR 
CLs and the optical [O iii] emission together with the similar kine- 
matics displayed by the optical low- and high-ionization lines 
jMazzalav eTai]|2O10l) suggest that the CLs and the low-ionization 
lines are produced by the same set of NLR clouds, or, else, that the 
NLR is clumpy on small scales (< 15 pc), with a mix of clouds at 
high and low density. If we assume that the former is the case, then 
we can set an upper value for the density of the CL gas as it must 
have densities lower than the critical density of the [O m] 5007 A, 
n^. = 7 X 10^ cm"^. That is, much lower than the critical density 
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of the CLs studied here (> 10^ cm"^^). This value is in complete 
agreement with the density of the CLR found to be between 10^ 
and 10* cm"^ by Rodriguez-Ardila et al. (in preparation), based on 
the analysis of optical [Fe vii] flux ratios measured for a sample of 
12 well-known AGNs. 

The emission of both low- and high-ionization lines (includ- 
ing the optical [O iii] 5007 A) are elongated along the direction of 
the jet observed in 6 cm with VLA, with an opening angle slightly 
larger than the one displayed by the jet (panel a of Fig. |4j. At 
the smaller scale structures traced by MERLIN (panels b and c of 
Fig. |4j the strongest NIR CL emission, particularly that of [Si vi] 
and [Si vii], coincides with the first portion of the jet, from SI to 
C. At point C, the jet bends from its original direction and re- 
gions of strong CL emission are now prominent at the NW and SE 
edges of the jet. This is compatible with the scenario suggested by 
iMacchetto et alj jl994h based on the [O m] images, where the jet 
cleans a channel trough the ISM, leaving behind hot ionized gas. 
As a result, the line emission is boosted along the edges of the jet 
where the gas is compressed. Additionally, cloud G is located in 
the region where the coUimated radio jet becomes lobe-like. 

Although the 6 cm radio emission is also observed towards 
the SW with VLA, the small-scale observations of MERLIN do not 
show evidence of emission in that direction. Pecontal et al. ( 1997) 
convolved the MERLIN data in order to obtain a spatial resolution 
comparable to that obtained with VLA, but did not detect emission. 
This suggests that the MERLIN observations are not deep enough 
to detect the SW jet structures. Therefore, it is not possible to make 
a small-scale comparison between the jet and the CL emission in 
these regions. Nevertheless, the spatial relation between the radio 
jet and the CL emission is remarka ble, suggesting tha t both are 
closely related. Indeed, very recently, IWang et ahl ( l20I2b presented 
exquisite Chandra X-ray images of the inner few hundred parsecs 
of NGC 1068 at a spatial resolution matching that of NIFS data. 
Extended X-ray emission is seen following that of [O m] HST, the 
CL emission (e.g. [Sivi]) and Merlin 6 cm. They found strong ev- 



idence of shocks, possibly related to the radio jet affecting knots 
B-l-C, cloud G and the region located 1.2 arcsec SW of the nucleus 
while the jet thrusts through clouds D-l-E and F without sings of 
strong interaction. This issue is explored further in Section |?!4l 

3.2.1 Comments on the morphology 

The H2 map (Fig. [3]l shows a rather different distribution from 
those of the other lines, displaying a ring-like structure around 
the nucleus, with the strongest emission located along the E- 
W direction at a distance of ~ 0.7 arcsec towards the NE. 
The same distribution is displayed by other H2 lines observed 
in NGC 1068, here we show the H2 2.12 jmi line as a exam- 
ple. The H2 emiss i on in NGC 1068 has been studied in detail by 
iGalliano & Alloid ( l2002h . These authors suggested that it traces 
a rotating molecular disc, consistent wit h the interferome t ric ob - 
servations of the '^C0(2-l) reported bv ISchinnerer et al.l fcOOOh . 
Additionally, iMulIer Sanchez et presented a study of 

the H2 emission of NGC 1068 from SINFONI^VLT integral field 
spectroscopy, focused on the central 1 arcsec. They interpret the 
data as evidence of hot gas, on scales of a few parsecs, fu- 
ell ing the AGN. The data p resen ted here are consisten t with that 
of iGalliano & AlloinI ( l2002h and lMuller Sanchez et ai] ( I2OO9I see 
their fig. 1), regar ding the size, shape and orienta tion of the H2 dis- 
tribution (see also Storchi -Bergmann et al ]|2012'). A more detailed 
analysis of the molecular gas using this same data set is carried out 
by Riffel et al. (in preparation). 

3.3 Kinematics of the CLs 

The CLs of NGC 1068 display very intricate profiles, with mul- 
tiple components and strong variations in their form and width 
from region to region. As an example. Fig. |5] shows the [Si vi] 
line profiles of spectra extracted from different spatial regions of 
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Figure 4. Flux distribution maps of [Si vi], in units of 10"'^ erg cm"^ sec"' . The contours correspond to the 6 cm emission observed with VLA (panel a) and 
MERLIN (pane l s b an d c) and [Oiii] emission (panels d and e). The labels SI, S2, C and NE indicated i n panel c conespon d to the radio knots defined by 
iGallimore et al.lfl996l) . The labels B— G indicated in panel e correspond to the emitting clouds identified bv lEvans et al.l fl99ll) . 



~ 0.2 X 0.2 arcsec^. Regions B-G were selected to approxi mately 
match the individual clouds identified bv lEvans et alj ( 1 199 ih in the 
[Oiii] image (see panel e of Fig. ID. 

In order to study the kinematics of the high-ionization gas sur- 
rounding the nucleus of NGC 1068, we have constructed channel 
maps of the strongest CLs (i.e. [S ix], [Sivi] and [Si vii]) observed 
in the NIR spectra. These maps were constructed by integrating 
the flux in small velocity bins (~ 150 km s"') along the line pro- 
files. Note that the high complexity displayed by the emission lines 
in this galaxy does not allow any meaningful kinematic analysis in 
terms of the standard velocity and velocity dispersion maps derived 
from the Gaussian fitting of the lines. 

Figs.[6]and|7]show the channel maps of [Si vi] and [Si vii], re- 
spectively. The central velocity of the integrated bin with respect to 
the systemic velocity of the galaxy is indicated in the upper right 
comer of each panel. The positive/negative values correspond to 
redshifted/blueshifted velocities. It can be seen that the channel 
maps of both lines are very similar. Some of the observed differ- 
ences can be attributed to the contamination of the [Si vii] line by 



sky line residuals and H2 2.48 fim emission (see section [STTT l. Bad 
sky subtraction is noticeable in the channels with Av = 154, -615 
and -1075 km s"', crossing the nucleus along the SE-NW direc- 
tion. Contamination by H2 is clearly distinguishable in the panels 
corresponding to Av = -1075 and -922 km s"', in the region lo- 
cated at ~ 1 arcsec E from the nucleus, where the H2 emission 
reaches its maximum (see Fig. |3}. Despite these differences, the 
agreement in the distributions is excellent. 

It is easy to see from Figs.|6]and|7]that both [Si vi] and [Si vn] 
are emitted by gas with a wide range of velocities, both positives 

and negatives, reaching values as high as 1500 and 1000 kms"'. 

The gas in the central regions (inner ~ 0.5 arcsec) is more turbu- 
lent, emitting at practically all the observed velocities. The line pro- 
files of these regions are broad, showing multiple peaks and strong 
wings. The CL emitting gas located further out from the nucleus 
is aligned in the NE-SW (PA~ 30°) direction, with a trend of the 
gas located towards the NE to show negative velocities while the 
gas located towards the SW displays positive velocities. This trend 
becomes stronger as we look further from the nucleus, with an in- 
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crease of the values of the observed velocities. Additionally, part of 
the gas located at the NE also shows positive velocities, so that the 
line profiles in these regions are split into two components, with 
one redshifted and the other one blueshifted (see Fig.|5}. This gen- 
eral pattern is not followed by the [Si vi] and [Si vii] emitting gas 
located at ~ 2.2 arcsec towards the SW from the nucleus, where 
the gas velocity is very close to that of the systemic velocity of the 
galaxy. Moreover, this system seems to be disconnected from the 
parcel of gas located next to it, towards the nucleus, because the 
latter displays much higher velocities. Overall, the [Si vi] channel 
maps are consistent with the ones presented by MS 11 (see their 
fig. 12), with the exception of the higher blueshifted velocity chan- 
nels, which are contaminated by H2 emission (see especially the 
two extremely strong blobs of emission in the -1400 km s"' ve- 
locity bin). 

The kinematics of the [Si vi] and [Si vii] emit ting gas derived 
in this work is compatible with that derived by iMazzalav et alj 
^2010) for the optical CLs along PA=22°. Although it is not pos- 
sible to make a point-to-point comparison between the two data 
sets (due to astrometric imprecisions), we conclude that both opti- 
cal and NIR CLs follow a similar behaviour: the nuclear lines show 
broad profiles with velocities close to the systemic velocity; fur- 
ther out, towards the SW, the strongest emission is produced by gas 
with receding velocities. Towards the NE, the line profiles are split 
into at least two components, with both receding and approaching 
velocities. 



Fig.[8]shows the channel map corresponding to the [S ix] emis- 
sion. It can be seen that this CL follows a similar trend as those 
of [Sivi] and [Sivn]. Note, however, that in the former the emis- 
sion appears more blurred than in the latter two due to a loss of 
spatial resolution in the y-band (see Section [3}. Furthermore, the 
[S ix] emission displays a smaller range of velocities than the sil- 
icon lines. The high-velocity gas (both approaching and receding 
components) in the silicon lines is mostly found in regions far from 
the nucleus (both towards the NE and SW), where weak or no [S ix] 
emission is observed. 

Overall, the integrated flux map of the CLs presents a clear 
asymmetry along the SW-NE direction (PA ~ 30°), that is, aligned 
with the radio emission observed with VLA (see panel a of Fig.|4ll. 
However, when we look at the channel velocity maps, it is possi- 
ble to discriminate two different behaviours: the high- velocity gas 
is extended along the SW-NE direction while the low-velocity gas 
(Av ~ -300-100 km s"') displays an asymmetry along the S- 
N direction. The latter corresponds to the original direction of the 
jet (formed by the radio knots SI and S2) before being deflected 
to its final direction (PA~ 30°). Assuming that the jet is initially 
aligned to the collimation axis of the ionizing radiation (i.e. the 
S-N direction) before bending to the NE, we expect that photoion- 
ization by the central source will be more relevant in that direction. 
In contrast, processes related to the interaction of the jet and its sur- 
rounding medium should be enhanced along the NE-SW direction. 
Taking this into account, it is possible to associate the low-velocity 
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Figure 6. [Si vi] channel maps. The centre of each velocity bin relative to the systemic velocity of the galaxy is indicated in the upper right comer of the 
corresponding panel. 



CL emission to gas photoionized by the AGN while gas radially 
accelerated due to the interaction with the radio jet should display 
the largest velocities. Indeed, Figs. [6] and [7] show that gas acceler- 
ated to velocities of up to 500 km s"' is present along the edges of 
the radio contours. This scenario would also account for the blue- 
and redshifted lines observed to the NE and SW of the nucleus, 
respectively. 

Currently, two main competing models have been proposed 
to explain the kinematics of the ionized gas in the nucleus of 
NGC 1068; on the one hand, the observed NLR kinematics are 
the resu lt of the interact i on of the nuclear gas and the radi o 
jet (e.g. lAxon et~Zlll99l ICapetti et alj [19971 : [Axon et alj|l998h . 
On the other hand, the observed kinematics can be explained by 
a biconical radial outflow, pointing to scenarios where the gas 



acceleration is produce by radiation and/or wind pressure (e.g. 
ICrenshaw & Kraemeil I2OO0I: | Das et al. Illooi). While these mod- 
els suggest very different scenarios, the observational data ob- 
tained to date are consistent wit h both models. Very recently, MS 1 1 
tested the scenario proposed bv lCrenshaw & Kraemeil ( I2OOOI) using 
[Si vi] channel maps derived from integral field spectroscopy SIN- 
FONl/VLT observations. Although the detailed description of the 
kinematic modelling method was left for a future publication, these 
authors find that the best-fitting model para meters are in agreement 
with those of previous long-slit HST studies jCrenshaw & KraemeJ 
l200d : lDasetaill2006h . Note, however, that the derived parameters 
may be aff^ected by the contamination of the H2 1.96 /im to the 
[Si vi] emission line discussed above, which particularly affects the 
blueshifted high-velocity channels. 
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Figure 7. Same as Fig|6] but for [Sivii]. Weak contamination by H2 emission can be seen in the panels corresponding to Av = -1075 and —922 km s"' 
at ~ 1 arcsec E from the nucleus. Additionally, bad sky subtraction afl'ects the channels maps in regions crossing the nucleus along the y-axis, especially 
noticeable at Av = 154 km s"'. 



The overall trend observed in the CLs channel maps derived 
from the NIFS observations is consistent with the main kinematic 
models proposed so far to explain the long-slit observations of the 
ionized gas in this galaxy. However, the richness of the 3D data pre- 
sented here gives a more global view of the nucleus of NGC 1068, 
and will guide future modelling. 

3.4 Emission-line ratios 

Emission-line ratios provide important information about physical 
properties of the emission-line gas (e.g. density, temperature and 
metallicity) averaged over all column densities along a particular 
line of sight. The excellent spatial resolution and broad wavelength 
coverage of the integral field spectroscopy data allow us to allevi- 
ate the limitations caused by projection effects by the analysis of 
the emission-line ratios produced by gas with different velocities. 
Although it is not possible to guarantee that parcels of gas with 
similar velocity are co-spatial, we can assert that parcels of gas 
with very different velocities are located in different regions along 
the line of sight. Thus, the resultant physical parameters are more 
representative of the different regions than the ones we would ob- 
tain integrating along the line of sight. In this section, we analyse 



the [Si vii]/[Si vi], [Feii]/Pa/3 and [Fen]/[Pn] emission-line ratios 
in order to gain information about the ionization structure and the 
presence of shocks in the inner regions of NGC 1068. 



3.4.1 The [ Si vii ]/[Si vi ] line ratio 

Few works in the literature have reported the simultaneous de- 
tecti on of [Sivi] and [S i vn] in NGC 1068. To our knowledge, 
onlv lMarconi et all ( 1 19961) and lRodriguez-ArdilT et al. (2006) mea- 
sured these two lines by means of long-slit spectroscopy. The rel- 
ative strengths of the silicon emission lines found in these two 
works pointed to a hard spectrum photoionizing source as the 
primary excitation for the majority of the emission. It was as- 
sumed that the source was th e nuclear central engine. Moreover, 
iRodrfguez-Ardila et al. I ( l2006h suggest that shocks generated by 
the radio jet NLR interaction could provide the additional power 
required for line formation. Individual [Si vi] and [Si vii] images, 
shown in Figs. [6] and [7| respectively, reveal resolved and extended 
emission in both lines, with at least three major separate compo- 
nents: one close to the nucleus, f ollowing the directi on of the radio 
jet; a second one at spot 'G' of lEvans et alj jl99ll) and a third at 
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Figure 9. [Si vii]/[Si vi] line ratio distribution measured for different velocity bins. The central velocity of each bin (relative to the systemic velocity) is 
indicated in the upper right comer of the corresponding panel. 



2 arcsec SW from the nucleus. Note that the latter two appear sig- 
nificantly stronger than or comparable to Bry emission. 

As both silicon lines are from the same element, their ratio is 
independent of the metal abundance. Because of their difference in 
ionization potential (AIP ~ 40 eV) the ratio can be used to map the 



ionization structure and mechanisms powering the CLs in greater 
detail than previous long-slit spectroscopy. 

Fig.[9]shows the [Si vii]/[Si vi] line flux ratio in NGC 1068 for 
the velocity bins between -1500 km s"' to 600 km s"' at a step of 
300 km s"'. Only ratio values with S/N>3 are plotted. It can be 
seen that a large variation in the value of the ratio is found across 
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the different regions, from ~ 0.4 to ~ 3, that is, nearly 1 dex. The 
largest values are found in an elongated region coincident with the 
inner portion of the jet before bending (S-N direction), where a ra- 
tio of ~ 1.9 is measured, and at ab out 1.3 arcsec SW from the centre 
(cloud HST-H: IWang et alj|2012!) . where a bright compact spot dis- 
playing a ratio of ~ 3 is observed. In velocity space, these two peaks 
correspond to gas at moderate blueshift (v = -300 km s"') and red- 
shift (v = 300 km s"'), respectively. The 'G' spot at ~ 130 pc NE 
from the nucleus appears as an isolated blueshifted feature with 
a [Si vii]/[Si vi] line flux ratio peaking at ~ 1.5. The smaller val- 
ues ([Si vii]/[Si vi] ~ 0.5) are predominantly found around the 'G' 
spot as well as in the emission region at 2.2 arcsec SW from the 
centre, and in some localized regions SW of the nucleus. In most 
other areas, the line r atio varies betwe en 0.9 and 1.4 . For reference, 
iMarconi et al.l 1 19961) and lRodriguez-Ardila et al.! ( i2006h measured 
a [Si vn]/[Si vi] integrated flux line ratio of 0.8 ± 0.3 and 0.6 ± 0.2, 
respectively. 

Taking into account uncertainties, the flux ratios measured in 
the GEMINI/NIFS data agree with the ones measured previously in 
the long-slit data except for the three peaks mentioned above. The 
fact that they do not coincide spatially with the nucleus, or are not 
symmetrically distributed around it, reflects the complexity of the 
inner regions of NGC 1068, where, for example, an anisotropic ra- 
diation field, gradients in density and/or selective absoiption could 
be s haping the observed ion ization structure. 

Iperguson et al.l ( 1 19971) present the results of a large number 
of photoionization simulations of the CLR in AGNs, with the ion- 
ization parameter U (//£] being the fundamental parameter of their 
models. They found that CLs fom at distances from just outside 
the broad-line region to ~ AQQL}^^^ pc, where L43 5 is the ionizing 
luminosity in units of 10"*^-^ ergs s"', in gas with ionization pa- 
rameter -2.0 < log U{H) < 0.75, suggesting that CLs will form 
close to the nucleus in high-density gas and further out from the 
nucleus in lower density gas. The peak equivalent width of each 
line is provided by the models. Since that quantity is referenced 
to the same point in the incident continuum, a comparison should 
indicate grossly the expected relative strengths of the CLs. Using 
table 1 of lFerguson et al] ( il997D . a [Si vii]/[Si vi] ratio of 1 .2 is pre- 
dicted. At first sight, this result agrees, within a factor of 2, with the 
values observed in NGC 1068. Note, however, that the calculation 
refers to the peak in the equivalent width distribution, which occurs 
at distances of a few parsecs from the central source and gas den- 
sity > 10* cm"^. At distances of about 100 pc, the ratio is expected 
to drop to a fraction of the peak value. 

The above result reinforces the need of an additional mech- 
anism powering the CL emission in NGC 1068. In the light of 
all observational evidence seen in the preceding sections, we are 
therefore let to consider the role of shocks, produced by interac- 
tions of either the radio jet or a radially accelerated outflow and 
the ISM. Fig. [9] shows, for instance, that spot 'G' is particularly 
enhanced in the blueshifted bins (v= -900 and -600 km s"'), 
with a [Sivii]/[Si vi] ratio of ~ 1.5 at the centre of the spot, 
surrounded by gas characterized by a line ratio < 1. As the 
NE region of NGC 1068 is very little affected by extinction (e.g. 
iKraemer & Crensha w 2000a,b), dust cannot be blamed for the in- 
crease of [Sivii] emission compared to that of [Sivi]. Moreover, 
spot 'G' is located at ~ 130 pc from the centre, which according 



The ratio of ionizing photon density to Hydrogen density, U{H) = 
0(H)/n{H)c, where <!>(//) is the flux of ionizing photons and c is the speed 
of light. 



to the results of lFerguson et al.l ( 1 1997b is far from the outer radius 
where [Sivii] emission is expected to occur (< 80 pc) assuming 
photoionization by the central source alone. Cloud 'HST-H' is also 
peculiar in this respect, as it shows the largest [Si vii]/[Si vi] ratio 
(~ 3) and is located at a similar distance as spot 'G' but in the 
opposite direction from the nucleus. Wang et al. (2012) found that 
both spot 'G' and cloud 'HST-H' also share similar properties at 
higher energies: both are conspicuous in X-ray emission and show 
significantly lower [Oiii]/X-ray ratio (~ 1) than the typical range 
spanned by photoionized clouds (~3-ll). They interpreted this as 
a higher collisional ionization at these locations and suggest that 
these clouds are strongly interacting with the jet, and that shocks 
driven into the obstructing clouds produce thermal X-ray emission. 
The increase in the [Si vn]/[Si vi] ratio at these two locations rein- 
forces this hypothesis. 

In this regard, it is very useful to check if models involving 
shocks are able to reproduce the ob served [Si vn]/[Si vi] ratios of 
NGC 1068. IContini & ViegasI l lioOll CVOl hereafter), for instance, 
presented a grid of simulations where the effect of shocks coupled 
to photoionization by the central source is employed to calculate 
the spectrum emitted by the NLR of AGNs. In this approach, the 
clouds are moving outwards from the galaxy centre, with the shock 
front forming on the outer edge of the cloud, whereas the ionizing 
radiation reaches the opposite edge that faces the active centre. The 
ionization due to the primary radiation (from the central source), 
to the diffuse radiation generated by the free-free and free-bound 
transitions of the shocked and photoionized gas, as well as the col- 
lisional ionization, are all accounted for The shock velocity Vs and 
the ionizing flux from the central source at the Lyman limit reach- 
ing the cloud, Fh (in units of cm"^ s"' eV"'), are the main input pa- 
rameters. Other parameters include the pre-shock density, iiq, and 
the pre-shock magnetic field, Bq. 

Tables 1-12 of CVOl show that shock-dominated clouds 
(F), = 0) with shock velocities in the range 300 - 500 km s"' 
strongly favour the production of [Si vii]/[Si vi] ratios larger than 
1.5 (models 56 and 70 for example). Moreover, model 70, with 
Vs = 500 km s"', produces a line flux ratio of ~ 3. For all these 
clouds, a solar metallicity, no = 300 cm"^ and Bo = lO"'' Gauss 
were adopted. When coupled to the presence of a strong radia- 
tion field from the central source (e.g. log Fh = 12, model 62), 
the value of the ratio drops to ~ 0.8. Thus, shocks locally enhance 
the emission of CLs even when taking into account its photoionized 
precursor As the maps of Fig. |9] show, blueshifted emission lines 
in the interval 300 - 900 km s"' characterize the regions where 
the largest [Si vii]/[Si vi] ratios are observed. That is, well within 
the range of shock velocities considered in the CVOl models. Note 
also that the gas is compressed downstream depending on Vs. For 
a Vs = 300 km s"' and no = 300 cm"\ the maximum downstream 
density n^ax is 1.04 x 10"* cm"^ while for Vs = 500 km s"', n^ax 
is 1.75 X 10'* cm"^ (see table 13 of CVOl). In both cases, n,nax is 
consistent with the upper limit for the CLR found in Section |T2| 

Additional support of the contribution of shocks to the 
CL emission in N GC 1068 can be obtained from the work of 
lAUen et al] ( 1200 8 h. They published an extensive library of fully 
radiative shock models that include both the radiative shock and 
its photoionized precursor. Similar to the models of CVOl, the 
shock velocity, the pre-shock density and the intensity of the mag- 
netic field are the main input parameters. Model predictions for 
the composite shock+precursor structure for the [Si vii]/[Si vi] ra- 
tio are found to vary from 0.5 to 1.2 for velocities between 500 
and 10 00 km s~', respectively , and no of 100 cm"^. Using this ap- 
proach, lExposito et al.l ( 1201 ll) were able to reproduce the CLs at 



© 0000 RAS, MNRAS 000, 000-000 



The CLR ofNGC 1068 13 



the different nodules in NGC 1068, which they attributed to a local 
jet-induced ionizing continuum. We should keep in mind, though, 
that the shock model predictions provide us with only a first-order 
approximation. In order to have a full description of the CLR based 
on model-fitting, we should suitably combine clouds under differ- 
ent physical conditions (i.e., clouds with different Vs and ng) with 
the proper weight so that the emission lines and the observed con- 
tinuum can be reproduced. This is far beyond the scope of this pa- 
per and is left for a future publication. What is most important here 
is that shock models are able to reproduce CL flux ratios at dis- 
tances as far as a few hundred parsecs, where photoionization by 
the central source predicts no or very faint CL emission. 

It is also very interesting to comment on the high-excitation 
channel with origin at the central engine and extending to the north, 
displaying a nearly constant [Si vii]/[Si vi] ratio of ~ 1.9 (see Fig.|9] 
at the V = -300 km s"' bin). This feature is coincident with the 
position angle of the radio jet before bending, continuing with that 
orientation beyond the point where the radio jet changes its ori- 
entation. At the distances from the centre at which this gas is lo- 
cated (£ I arcsec), photoionization by the central source should be 
enough to power the observed CLs. Not e, however, that fa int X- 
ray emission is also seen in the images of lWang et alj ( 1201 2h along 
the high-excitation channel, revealing the presence of hot plasma. 
Shocks along this region should not be discarded. 

All the observational evidence above points out that 
shocks must contribute significantly to CL emission, confirm- 
ing early suggestions base d on long-slit observations (e.g. 
iRodrfguez-Ardila et al.l |2006|) . Previous works had considered 
this possibility but overall found that photoionization by the 
central source is the main driving m e chanism of the CL 
emission (e.g. 'Oliva & Moorwoo 3 '1990'; 'Marco ni et J] 1 19961 : 
aemer & Cr enshaw 2000b; Mazzalay et al. 2010). The evidence 
shown here observationally confirms the relevance of shocks as an 
important component of the CL emission in NGC 1068 and, thanks 
to the spatial and spectral coverage of our integral field spectro- 
scopic data, we are able to identify the particular regions where 
shocks play a significant role. A comparison with the exquisite X- 
ray map of the inner 200 pc of that galaxy is in accord with this 
scenario. 

In summary, the [Sivii]/[Si vi] line flux ratio map found for 
NGC 1068 reveals localized regions where the gas excitation can- 
not be easily explained by simple photoionization by the cen- 
tral source. Moreover, the fact that these regions follow closely 
the radio morphology at some specific locations, invokes a sce- 
nario where shocks originating from the interaction of the emit- 
ting clouds with the radio jet contribute to enhance the ionization 
state of the gas. However, other sources of shocks should also be 
considered (the interaction of outflowing clouds with the ISM, for 
instance). Models taking into account the effect of shocks, with ve- 
locities similar to those measured in NGC 1068, can account for the 
high [Si vii]/[Si vi] line flux ratio observed at some locations in the 
NLR of that object. 

3.4.2 The role of shocks in the formation of [Fen]: the 
[Fe II ]/Pa/3 and [Fe n ]/[P ii ] line ratios 

The origin of the [ Fen] emission in galaxies has been ex ten- 
sively studied (e.g. iForbe s & Ward 1993; Simpson & Mea dows! 
1998; Rodn'gu ez-Ardila et al.ll2004 120051; IStorc hi-Bergm ann etai] 
2009)- However, the excitation mechanisms and location of the 
[Fen] emitting gas are still a subject of debate. Due to the low- 
ionization potential (16.2 eV), the general consensus is that [Feii] 



is formed in a partially ionized zone. Photoionization by X-rays 
and shocks have been suggested as the most important mechanisms 
contributing to the fomation of the [Fe ii] in this region. 

In this context, the [Feii]/Pa/3 ratio is considered as an im- 
portant tracer of the origin of the [Fe n] emission. There is an in- 
creasing progressi on in this ratio from pure p hotoionization to pure 
shock excitation jSimpson & Meadowsl [l998) as one goes from 
starbursts (displaying typical values of [Fen]/Pa/? ~ 0.5) through 
AGNs to SNRs and shocks (with values higher than ~ 2). There- 
fore, this line ratio can be used to measure the relative contribution 
of photoionization and shocks in the formation of the [Feii] lines. 
Fig. 1 101 shows the [Fen]/Pa/? line ratio calculated for four different 
velocity ranges in the nucleus of NGC 1068, with the 6 cm VLA 
contours overlaid. Low values of [Fe iiJ/Pa/? (< 2) and a homoge- 
neous distribution characterize the gas with velocities in the range 
Av = -300 ± 150 km s"' . On the other hand, the gas with velocities 
close to the systemic velocity and higher show a very different line 
ratio moiphology, resembling partially concentric rings with a sym- 
metry axis coincident with the radio jet. The values of [Fe ii]/Pa/3 
increase in regions further out from the nucleus, reaching values in 
the SNRs range (> 2) in the regions associated with the NE end of 
the jet. Both the shape and values of the [Fen]/Pa/3 strongly suggest 
the presence of shocks in the inner regions of NGC 1068. These 
contribute to the fomation of [Fe n] in the outer regions while pho- 
toionization probably dominates closer to the nucleus. 

Another useful tool to constrain the or igin of [Fe n ] lines is the 
[Fell] 1.26 fJ.mJ[Fa] 1. 19 urn line ratio. lOliva et all ( 1200 1[) pro- 
posed that in objects with low [Fen]/[Pii] ratios (< 2), shocks do 
not play an important role in the line excitation while large values 
of this ratio (> 20) indicate that the emitting gas has recently passed 
through a fast shock and it is likely that the lines are produced 
by shock excitation. Although the [Pn] emission in NGC 1068 is 
weak, it was possible to construct [Fe ii]/[P n] maps for the inner- 
most regions. These maps are presented in Fig. [TT] The measured 
values of this line ratio vary between ~ 0.5 - 7.8 and, as in the case 
of the [Fen]/Pa/? ratio, its distribution shows a spatial correlation 
with the radio jet. Alth ough the measured values are lower than the 
value of 20 derived bv lOliva et al.l ( 1200 ih for shocked gas, the re- 
gions associated with the NE component of the jet are much higher 
than the lower limit of 2 suggested by these authors, indicating an 
important increase in the relative abundance of Fe/P due to shocks. 

Additional evidence for both photoionization and radio 
jet/shock excitation of [Fen] in the nuclei of AGNs can be found, 
for example, in the work presented by Storchi-Bergma nn et alj 
(2009, 2010). These authors analysed NIFS/Gemini integral field 
spectroscopy data of the well-known Seyfert 1 galaxy NGC 4151, 
finding a strong influence of the radio jet present in the centre of 
this galaxy on the [Fen] emission. From a similar analysis to the 
one presented here, these authors concluded that photoionization 
dominates in the inner region while shocks dominate in the outer 
region of the galaxy. 

A more detailed discussion on this subject, together with the 
analysis of the low-ionization lines of NGC 1068, is carried out by 
Riffel et al. (in preparation) using the same set of data. 

In summary, our findings indicate that shocks produced by the 
interaction between the radio jet and the ambient gas are respon- 
sible for most of the [Feii] emission observed in the nucleus of 
NGC 1068, especially in the regions located at the NE of the nu- 
cleus. 
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Figure 10. [Fe ii]/Pay6 line ratio distribution measured for different velocities bins. The central velocity of each bin (relative to the systemic velocity) is indicated 
in the upper right comer of the corresponding panel. 
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Figure 11. [Feii]/[Pii] line ratio distribution measured for different velocities bins. The central velocity of each bin (relative to the systemic velocity) is 
indicated in the upper right comer of the corresponding panel. 



4 SUMMARY AND CONCLUSIONS 

In this paper we have presented new J- and A^-band 2D spectra of 
the inner ~ 300 pc of NGC 1068. The data were obtained with the 
Gemini NIFS integral field unit spectrometer, which provided us 
with high spatial and spectral resolution sampling. 
Our main results can be summarized as follows. 

(i) The CLs observed in the spectra of NGC 1068, [Sivi], 
[Sivii], [Six], [Caviii] and [Alix], display very asymmetric and 
inhomogeneous distributions. The emission is elongated along the 
NE-SW direction, with the stronger emission preferentially local- 
ized to the NE of the nucleus. We observe CL emission up to a 
distance of ~ 170 pc to the NE and SW of the nucleus. 

(ii) Hydrogen and [O iii] low-ionization lines share a similar dis- 
tribution (bright knots and arcs) as the CLs, suggesting a common 
origin. On the other hand, the morphology displayed by the molec- 
ular hydrogen lines is completely different. The [Fe ii] emission ap- 
pears to have a shock-excited component. 

(iii) CLs are emitted by gas covering a wide range of velocities, 
up to ~ 1000 km s"' relative to the systemic velocity of the galaxy. 
There is a trend for the gas located on the NE side of the nucleus to 
be blueshifted while the gas located towards the SW is redshifted. 
Nevertheless, some redshifted emission is observed NE from the 
nucleus, consistent with the double-peak profiles observed in opti- 
cal high- and low-ionization lines. 

(iv) The comparison of the line flux distributions with the ra- 
dio emission suggests that the morphology of the emission-line 
gas is closely related to the radio jet observed in the nucleus of 
NGC 1068. On the other hand, although the overall gas kinematics 



is consistent with being driven by the radio jet, other mechanisms 
are also possible. 

(v) The ionization structure, mapped by the [Si vii]/[Si vi] 
emission-line ratio, is highly intricate. While this ionization map is 
not easily explained by photoionization by a central source, models 
taking into account the effects of shocks are able to reproduce the 
high values of [Si vii]/[Si vi] observed in some regions. 

(vi) The presence of shocks in the nucleus of NGC 1068 is sup- 
ported by the high values and distribution found for the [Fe n]/Pa/8 
and [Fe n]/[P ii] line ratios. 

The richness and complexity of the structures we observe, 
both in terms of velocity fields and emission-line ratios are remark- 
able. Therefore, the use of integral field spectroscopy is essential 
in order to map the circumnuclear environment of nearby AGN, to 
disentangle the effects of photoionization and shocks on emission- 
line properties, and to study the effects of radio jets on their gaseous 
environment. 
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